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SUM1ARY 


This progress report summarizes the research work performed at the 
Catholic University of America on the research grant entitled 'Distributed 
Active Control of Large Flexible Space Structures,” funded by NASA/Goddard 
Space Flight Center, under the grant number NAG 5-749, during the period of 
March 15, 1986 to September 15, 1986. 

In this report we first update the research work relevant to the 
project. Then the research accomplished during the above period will be 
stated. The report is then concluded by a discussion of current and future 


research work. 
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Introduction 


The advent of a space transportation system, such as the space shuttle 
makes it possible to conceive of very large satellites and spacecraft which 
could be carried into space and deployed, assembled, or constructed there 
for such diverse purposes as communications, surveillance, astronomy, space 
exploration, and electric power generation [2]. These large space 
structures (LSS) concepts range from central rigid bodies, to the solar 
electric propulsion spacecraft and the generic polar platform. Two control 
problems for LSS are attitude control and shape control [24]. The former 
involves maintaining a given orientation of the spacecraft, e.g., with 
respect to the sun or earth; the latter involves maintaining the shape of 
critical structures of LSS. LSS are distributed parameter systems that 
possess many low resonant frequencies and have very stringent requirements 
for shape, orientation, alignment, vibration suppression and pointing 
accuracy. These requirements lead the control designers to the concept of 
active control of LSS with various sensors and actuators located about the 
structure and operating through on-line computer controllers to tailor the 
performance and behavior of the system. 

There has been considerable interest in the area of active control of 
LSS [l]-[24], A number of control schemes were proposed for large flexible 
space structure (LFSS), but they all represent one form or another of modal 
control [11]. Two main modal control schemes are the coupled Modal Control 
and Independent Modal Space Control (IMSC). The former employs an active 
controller consisting of a state estimator and a state feedback law; the 
latter controls each mode independently by means of the modal filter [13]. 
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In this report we first update the research in active control of LFSS 
to provide a mathematical framework for the project. Then the modelling of 
a generic polar platform is developed using finite element method. After 
that some project research results are presented. We then discuss the 
current and future research effort. A conclusion will summarize the report. 
2. Updated Research 

In this section we first present the mathematical description of large 
space structures (LSS) and then discuss the two main control schemes for 
this type of structures. 

2.1 Mathematical Description of LSS 

The LSS may be described as a continuum by the following partial 
differential equations 

2 

M(P) 9 U 2 P,t '- + L u(P,t) = f(P,t) (2.1) 

at 

where u(P,t) * displacement of an arbitrary point P 

L * Linear differential self-adjoint operator of 

order 2p, expressing the system stiffness. 

M(P) * distributed mass 
f(P,t) = distributed control force 

Equation (2.1) must be satisfied at every point P in the domain D. The 
displacement of u(P,t) is subject to the boundary conditions 

u(P,t) = 0 for i=l,2,...,p (2.2) 

where B. are linear differential operator of order ranging from 
0 to <2p-l). 


5 



i 


The associated eigenvalue problem is formulated by 

I* O r (P) = X r M(P) * r (P) (2.3) 

for r=l,2, . . . ,o o 


with the boundary conditions 

B i $ r (P) = 0 (2.4) 

for i=l,2,...,p; r=l ,2 , . . . , °0 

where x^is the rth eigenvalues and <t>^(p) is the eigenfunction associated 

with x . Sometimes $ (P)is called the mode shape. 

r r 

Equations (2.3) and (2.4) can be solved to obtain the solutions of \ 
, and in addition, if the operator L is positive definite, then all 
eigenvalues are positive. The mode frequency (or natural frequency) is 
defined as 

i 

u> r = ( X r ) ; r=l,2, ... (2.5) 

using expansion theorem [11], the solution of u(P,t) can be obtained as 

oo 

u(P,t) = 4> r (P) U r (t) 

r=l 

where u r (t) satisfies 


- 6 - 


and 


( 2 . 6 ) 



u (t) + u> u (t) = f (t) 

r r r r 


-r (t> ■ j r 


(P) f(P,t) dD 


( 2 . 8 ) 


In practice, the infinite series in (2.6) is truncated as 


u(p,t) = y u (t> 

r r 

r=l 


(2.9) 


where M is chosen to be sufficiently large so that u(P,t) can be represented 
with good fidelity. 

Since M may be quite large, it is not reasonable to control all M 
modes. Hence we select N modes to control. These N modes are called 
Controlled Modes . The remaining R modes (R 3 M-N) are called Residual Modes 
(uncontrolled modes). 

Equation (2.7) can be transformed into the state equation form as 
follows 


*c (t) = A c x c (t) + f c (t) 

= A R X R (t) + f R (t) 
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( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 
(2.17) 


7 


PRECEDING PAGE BLANK NOT FILMED 


Discrete Actuators ; Since it is impossible to implement distributed 


control forces, the distributed control force f(P,t) is realized by k 
discrete point force actuators 

k 

f(P,t) = Y 5 (P-P . ) F.(t) (2.18) 

i=l 1 1 

where 6(P-Pi) is a spatial Dirac Delta function. 

Substituting (2.18) into (2.8) yields 




k 






f r (t) 

- I 
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r l 

At) ; r 

=1,2,. 

(2.19) 

Now 

using (2.19) we obtain 





f c (t) = B c F(t) 




(2.20) 

and 

F R (t) = B R F(t) 




(2.21) 




w ••• 

w 





*2 (P 1> 

W ••• 

VV 




B c 

* 

* 

• 


(2.22) 



vv 

V P 2> ... 

VV 






°Nxk 









1 




< * , N+l (P l ) 

$ N+1 (P 2 ) 

*N+1 ( V 




^N+2 (P 1 ) 

WV 

' * ’ WV 



b r = 

. 

• 


• 

(2.23) 



VV 

V P 2> 

. . . $ 

m'V 





°Mxk 







8 


and F(t) = [Fj Fg .... F R ]' ( 2 . 24 ) 

Discrete Sensors : Suppose there are s discrete sensors consisting of p 

velocity sensors, located at p locations P yl P^ 2 , . .., P yp and q 
displacement sensors, located at q locations P^, P^, ..., P^. Using 
(2.9) we can express the velocity sensor output as 

M 

Yn = = £ W "r (t) 


V]' 


r=l 


( 2 . 25 ) 


for j = 1,2,..., p, and the displacement sensor output as 

Vi = U(P Dj' t) = V P Dj> U r <t > 

for j = 1,2,..., q. 

If we define the sensor output vector as 


y (t) = i*i y 2 ••• y p Vi y p+2 

then using (2.25) and (2.26), we can write 

y(t) * C Q x c (t) + C R x R (t) 


Vq' 


where 


C C * bdIC CV C CD ] 
C R 8 bdlC RV C RD ] 


and 
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( 2 . 26 ) 


( 2 . 27 ) 


(2.28) 

(2.29) 

(2.30) 


( 2 . 31 ) 


( 2 . 32 ) 


( 2 . 33 ) 
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The development in this section can be now summarized by the following modal 
state equation in matrix form 
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- 



1 




— 

A c 

0 



+ 

B c 

*R 


0 

a r 


X R 
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_ 
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F ( t) (2.34) 


Y (t) = [ C C C R J 


(2.35) 
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2.2 The Coupled Modal Control Scheme 


In this scheme, the active controller consists of two parts: 

a) a state estimator that accepts the sensor output y(t) and produces 
an estimate Xg(t) of the controlled state Xg(t) and 

b) a linear state feedback control law that multiplies the estimated 
controlled state Xg(t) by constant gain to produce the actuator input F(t). 

The state estimator can be of the Luenberger type that is described by 

CIS] 


x c (t) » <A C - G c C c )I c (t) + B c F(t) + G c y(t) (2.36) 

The state estimator error is defined as 
e c (t) = 5 c (t) - x c (t) (2.37) 

Therefore 

e c (t) = xit) - i c (t) = (A c - G c C c )e c (t) + G c C R x R (t) (2.38) 

The linear state feedback law is defined by 

F(t) = K c x c (t) (2.39) 

Now substituting (2.39) into (2.34) and (2.38), we obtain the composite 
state equation as follows: 


*c (t) 


V B c K c B c K c 0 
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e c (t) 


0 Wc G cS 


(D 
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ft 

> <t} 


b r k c Vc a r 

k ^ 


> (t, j 


The objective of the couple modal control scheme is to select Kg and 6g such 
that the composite system (2.40) is asymptomatically stable. In (2.40) the 
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control spillover is represented by and the observation spillover is 
represented by C R . The coupled modal control scheme is illustrated in 
Figure 1. 

2.3 The Independent Modal Space Control Scheme (IMSC) 

The IMSC scheme was developed by Meirovitch and others [10] to control 

each mode independently. First an auxiliary variable V^t) is defined by 

v (t) = 1 /oj u (t) (2.41) 

r r r 

Therefore using (2.7) we can write 

v (t) =1/(0 u (t) = - U) 2 u (t) + f (t) (2.42) 

r r r r r r 

Now using (2.41) and (2.42), we obtain 


x r (t) = A r x r (t) + W r 


where 


A 

r 





w r = l 0 


T 

V“rl 


and x r = [u^ 




T 


(2.43) 

(2.44) 


(2.45) 

(2.46) 


The essence of the IMSC is to choose W r such that W r depends only on 
alone: 


W = W ( oi ) (2.47) 

r r r 

The optimal design of W^, was discussed in [12]. 


Modal Filter : To implement the IMSC scheme, the modal displacements 

u f (t) and modal velocities u r (t) are required. A device, called modal 
filter accepts the measurements of displacement u(P,t) and velocity u(P,t) 

at every point P at all time and produces the modal displacement u r (t) and 
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modal velocity u r (t). The mathematical representation of a modal filter is 
given by 


u r (t) = I M(P) <t r (P) u(P,t) dD (2.48) 

u (t) = [ M(P) 4> r (P) u(P,t) dD (2.49) 

r D 

It is noted that distributed sensors are needed to obtain the 
displacements u(P,t) and velocities u(P,t). 

Implementation of IMSC using discrete activators and sensors 

The implementation of IMSC requires distributed control, i.e. Control 
forces are applied at every point of the structure, and also requires 
distributed sensors, i.e. the displacements and velocities are measured at 
every point. This is however not realizable in practice. Thus discrete 
actuators and discrete sensors are used instead. 

Let us consider the problem of controlling N modes by means of discrete 
actuators. The actuator forces can be treated and distributed by writing 

N 


f (P,t) 



6 (P-Pj) 


(2.50) 


substituting (2.50) into (2.8) yields 

N 

f (t) = E* (pj) f (t) 

r -;=i r J 


(2.51) 


Introducing the matrix B as 


B = 


(P 1 > * 1 (P 2 ) 

W VV 


VV V p 2> 


w 

VV 


W j 


(2.52) 


we obtain 


f(t) = B F(t) 


(2.53) 
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where 


(2.54) 


f(t) = (f 1 f 2 .... f N ] T 
F(t) = [F 1 F 2 .... F n ) T (2.55) 

Suppose B is nonsingular, then the actual control vector F(t) can be 
synthesized from the generalized control vector f(t) by 

F(t) = B -1 f(t) (2.56) 

Next let us assume there are m sensors capable of measuring 
displacements and velocities at the discrete points P=P k (k*l , 2 , . . . ,m) . 

Then from these measurements the entire displacement pattern u(P,t) and 
velocity pattern u(P,t) can be estimated by using various interpolation 
functions. The structure of the IMSC scheme using discrete actuators and 
sensors is shown in Figure 2. 

2.4 Evaluation of the control schemes 

The design of the coupled modal control scheme is straight forward 
since it is simply a pole placement problem if the observation spillover is 
negligible. Indeed if we set C R =0 in (2.4), then the eigenvalues of the 
composite system are the eigenvalues of (Aq+ b c k^ 5 (A c ~G c C c )and A R due to 
block triangularity. There are several techniques such as in [15,16 ] that 
help the control designer to select KC and G^ such that A^+B^K^ and A^-G^X^ 
have desired eigenvalues. The control spillover is represented by 
x R (t) = A R x R (t) + B R K c [x c (t) + e c (t)] (2.57) 

Since e^(t) and x^(t) will decay to zero, the control spillover will cause 
unwanted excitations of the residual modes, but it cannot shift the residual 
mode frequencies. Consequently the control spillover results in some 
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unwanted oscillations of the residual modes. So when the observation 
spillover is absent, then the control spillover degrades the system 
responses but cannot destabilize the system. 

When observation spillover exists, then instability may occur when some 
of the eigenvalues of the composite system (2.40) are unstable. 

The control and obsesrvation spillover can be minimized if the 
actuators and sensors are located at (or very near to) the zero of the mode 
shapes of the residual modes. This concept is limited since it creates 
uncontrollable and unobservable triple (Ag, Bg, Cg) that hinders the design 
of state estimators and state feedback laws. In addition, the freedom to 
locate sensors and actuators arbitrarily is rarely available to the designer 
since these locations are often already determined by structural 
considerations. 

It is well known [12] that the IMSC scheme is capable of eliminating 
control and observation spillover. However it requires the implementation 
of distributed actuators and distributed sensors that unfortunately cannot 
be provided by the current state-of-the-art. Consequently, discrete 
actuators and discrete sensors have to be used. A question that is still 
not clearly answered is how many discrete actuators and sensors are 
necessary for a successful implementation of the IMSC scheme. In addition, 
the IMSC scheme when being implemented by discrete actuators also has 
control spillover into the uncontrolled modes. However as pointed out in 
[12], the control spillover effect is not very important. 

Both control schemes, the coupled modal control and the IMSC require 
that a closed-form solution of the eigenvalue problem exist. Unfortunately 
the vast majority of continuous systems leads to eigenvalue problems that do 
not have closed-form solutions, owing to nonuniform mass or stiffness 
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distributions. Hence in this case, it is necessary to seek approximate 
solutions of the eigenvalue problem. 

3. Modeling of the Polar Platform 

One of the project activities is to obtain a model for a generic polar 
platform whose physical parameters are given in Table 1. 

A NASTRAN finite element analysis package (McNeal Schwendler version) 
was employed to model the platform and its appendages by 52 modes. The 
model is illustrated in Figure 3. The platform itself was modeled by 28 
nodes, the solar panel and the astromast by 16 nodes and the engineering 
module by 8 nodes. 

In the developed model, the payload carried by each module was assumed 
to be distributed uniformly throughout all nodes. Also the engineering 
module was assumed to be very rigid as compared to the platform or the solar 
panel. 

The model considers only those degrees of freedom of the platform that 
contribute to its bending and torsional vibrations. 

Table 2 lists the first ten modes of vibrations. The first two natural 
frequencies are the first rigid body modes in pitching and rolling. The 
remaining modes are the flexible modes of the platform. 

It can be seen that the flexible modes are very low and highly 
condensed over a narrow band between 0.054Hz and 0.448 Hz. A sample of the 
first three modes of vibrations are shown in Fig. 4-6. 

4. Project Research Results 

In this section we will present some research results obtained by the 
principal investigators during the research period stated in the summary. 

One of the project objectives is to evaluate several candidate control 
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schemes in order to select the most appropriate control scheme for the LFSS. 
The study of coupled-modal control and Independent Modal Space Control 
involves in such problems as pole allocation, eigenvalue assignment, state 

estimator design, canonical transformation etc During the 

investigation of the above topics, some results have been found and are 
summarized below. 

4.1 Canonical Transformation 

Canonical transformation has been proven to be very useful in the 
design of state estimators and state feedback. In [14] two canonical 
transformations were developed for a class of time-varying multivariable 
control systems. Since most physical systems such as LFSS are time-varying 
in nature, this type of time-varying canonical transformations will be very 
helpful in design of more practical control schemes for LFSS. 

4.2 Eigenvalue Assignment and Pole Allocation 

The eigenvalue assignment or pole allocation problem is very essential 
to the control of LFSS since it is related to the system stability. Two 
common problems in control system design for LFSS are the selection of a 
state feedback gain to shift the natural modes of the LFSS to a set of 
desired damping modes, and the design of a state estimator whose eigenvalues 
can be arbitrarily assigned. These two problems can be treated as 
eigenvalue assignment or pole allocation problems. In [16] it was shown 
that eigenvalues can be arbitrarily assigned to a class of time-varying 
multivariable systems. Canonical transformation developed in [14] was 
employed to design the state feedback gain. The distributed control problem 
was considered in [17] where a control scheme was developed to arbitrarily 
allocate closed-loop poles to a distributed time-invariant system that is 
controlled by several pairs of sensors and actuators. 
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4.3 State Estimator Design 


State estimator design is very crucial to the control of LFSS since the 
controlled state vector has to be estimated for the coupled modal control 
case and the modal state for the IMSC case. In [18] a new algorithm was 
proposed to design reduced-order state estimators for a class of 
time-varying multivariable control systems that are uniformly observable. 

The result in [18] is an extension of that in [15] where a full-order state 
estimator was designed. 

4.4 Simulation of Control Schemes 

Since the generic polar platform is relatively a complicated structure, 
we start the study of the control scheme by considering the control of a 
simply supported beam with one actuator and one sensor. Complete 
understanding of the control of this simply supported beam will help us to 
gain some insight of the control of the complete structure. The beam 
dynamics are modelled by the Euler-Bernoulli partial differential equation 

2 

3 34 

m — -- u(x,t) + El 2 u ( x * fc ) = f(x,t) (4.1) 

at ax 

For simplicity we set the mass m, the moment of inertia I, the modulus of 
elasticity E and the length of the beam to unity. The boundary conditions 
for this simply supported beam are 

u(0,t) = u(l,t) = 0 (4.2) 

a 2 a 2 

---u(0,t) = ---U(l,t) = 0 (4.3) 

ax ax 
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(4.4) 


The eigenvalue solutions for this case are 

2 

u> k = (kit) 

= sin (kirx) 


(4.5) 


The beam is controlled by a single point actuator at x = 1/6: 

f(x,t) = (1/2) * 5 6 (x - 1/6) f(t) (4.6) 

and the displacement is measured by a single point sensor located at x = 
5/6: 

y(t) * u(5/6, t) (4.7) 


Using the Matlab and the control tool box software package, an active 
controller was designed to control the first three modes of the beam. To 


minimize the following unweighted energy 

.eO 


J » J [E( t) + 0.1f fc (t)]dt 

o 

E(t) = 1/2 


where 


3 

E 

k=l 


, 12 . 2 
( “ k \ + \ > 


the control gain was determined to be 

K c » (-3.87 -3.86 -3.86 -1.01 -6.08 -20.3] 


(4.8) 


(4.9) 


(4.10) 


A state estimator was designed to estimate the controlled modes. The 
state estimator gain Gg given below: 

T 



will assign stable eigenvalues to the state estimator. Several simulation 
runs were made to study the performance of the active controllers and the 
effect of the control and observation spillover. 

In all simulation runs, an impulse disturbance was created to excite 
the system. Figures 8 and 9 present the effects of control and observation 
spillover. Me note that when spillover exists, the displacement did not 
reach the zero steady state. It has an oscillation with small amplitude 
caused by control spillover into the uncontrolled modes. Even though the 
observation exists in this case, the system is till stable. Fig. 10 and 11 
display the sensor output and actuator input when no spillover exists. Me 
notice that the displacement was decayed to zero. The control spillover 
into the 4th and 5th uncontrolled modes is illustrated in Fig. 12 and 13, 
respectively. The eigenvalues of the composite system [Equation (2.40)] for 
all possible cases of spillover are given in Table 2a. Table 2b shows the 
eigenvalues of the closed-loop system, the state estimator, controlled 
subsystem and uncontrolled subsystem. 

In Table 3a we note that when no observation spillover exists, the 
eigenvalues of the composite system consist of the eigenvalues of the 
closed-loop system, the state estimator, and the uncontrolled subsystem. 
This agrees with the discussion made in Section 2.3. 

5. Current and Future Research 

Continuing the research stated in Section 4, we currently focus on the 
study of spillover minimization. As we have seen, the system simulation 
results showed that the system performance is degraded at the presence of 
control spillover. The LFSS can also become unstable if the observation 
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spillover exists. Therefore spillover minimization study is very crucial to 
the control of LFSS. Me are investigates the minimization of the spillover 
by 

(a) Proper locating of actuator and sensors. 

(b) Determining a sufficient number of actuators and sensors needed. 

(c) Applying the IMSC method to the generic polar platform. 

The above study is carried out mostly by simulating the system for 
various situations using system simulation languages in order to set up a 
guideline for selecting and locating actuators and sensors. 

After the control of the simply-supported beam has been studied 
throughout, we will apply the resulting control scheme to the complete 
generic polar platform. Several control issues such as given below should 
be addressed in the future research: 

(a) Evaluation of the application of coupled-modal control and IMSP to 
the generic polar platform in terms of feasibility, reliability, spillover, 
and computational effort. 

(b) Implementation of an adaptive control scheme for the system. 

(c) Study of a robust controller when the system parameters are not 
well-known and when unpredictable disturbances are expected. 

6. Conclusion 

In this report we first updated the research related to project. The 
modelling of the generic polar platform using finite element method was 
discussed. Me then presented some preliminary results of the control system 
study using computer simulation and some theoretical research results of 
such problems as eigenvalues, state estimator, and state feedback design. 
Current research effort was discussed and future research activities were 
assessed. 
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FIGURE Is The Structure of the Coupled Modal Control 
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FIGURE 3 


CO 



: The Finite Element Model of the Generic Polar Platform 
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POLAR PLATFORM- 



FIGURE 4: The First Mode of Vibration 
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HOSE 1(0 HZ) 
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HOPE 3 (9.0515 HZ) 
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FIGURE 10: Actuator Input with no Control and Observation Spillover 
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VALUE 


PARAMETER 

OVERALL LENGTH OF THE PLATFORM 

PLATFORM WIDTH 

PLATFORM HEIGHT 

LENGTH OF SOLAR ARRAY 

WIDTH OF SOLAR ARRAY 

NUMBER OF MODULES 

MASS/MODULE 

LENGTH OF ENGINEERING MODULE 
HEIGTH OF ENGINEERING MODULE 
MASS OF ENGINEERING MODULE 
LENGTH OF STRUCTURAL MEMBERS 
STRUCTURAL MEMBER OUTER DIAMETER 
STRUCTURAL MEMBER INNER DIAMETER 
STRUCTURAL MEMBER MATERIAL 
DESITY OF THE ASTROMAST SUPPORTING 
THE SOLAR ARRAY 
ASTROMAST FLEXURAL RIGIDITY 
ASTROMAST TORSIONAL RIGIDITY 


45 ft 
6 ft 
6 ft 
60 ft 
15 ft 
6 

600.00 kg 
9 ft 
6 ft 

5000.00 kg 
6 ft 

1 in. 

0.875 in. 

Aluminium 

2.2 kg/ft 
15.20 xlO^lb in 
& 3u-H^xl0 6 lb in 


TABLE 1: Physical Parameters of the Generic Polar Platform 
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1— 
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io" 1 
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io -1 
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TABLE 2: The First Ten Modes of Vibrations 
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TABLE 3a: Eigenvalues of the Composite System 
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